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Abstract 


The delivery of inventory to customers involves the use of vehicles for which purchase or lease agreements may 
need to be signed months or even years before the start of actual delivery operations. Such long lead times make it 
impractical to build a fleet incrementally and motivated study of the Strategic Inventory Routing Problem (SIRP). 
SIRP focuses on estimating, in advance of the start of actual delivery operations, the minimum size (or cost) vehicle 
fleet required to supply inventory from a central depot to spatially dispersed customers. A fundamental! difference 
between SIRP and the associated tactical IRP’s (route an existing fleet, given knowledge of actual customer 
inventory levels) is that all possible realizations of the tactical problem must be considered, at least implicitly, when 
solving SIRP. This paper generalizes Larson’s approach for SIRP through the use of period and phase of customer 
replenishment as additional decision variables. Routing solutions based on customer-specific period and phase of 
replenishment are developed for a simple model of the tactical routing problems the fleet will eventually encounter. 
Estimates of the fleet size required are developed on the basis of these routing solutions. The period / phase 
approach can be generalized to take long-term operating costs into account. Computational tests show that the new 
approach yields significant reductions in solution cost when the vehicle is large enough to replenish several 
customers in a single trip and/or when there is significant variation in the maximum inter-replenishment intervals of 
the customers. 


Keywords: Strategic Inventory Routing Problem 


Inventory Routing Problems (IRP’s) arise when both inventory and routing considerations are 
incorporated in a model to adequately capture the characteristics of a distribution system. Typical 
routing considerations include constraints on route length and/or duration, vehicle carrying capacity, 
vehicle / location compatibility and the times at which deliveries may be made. Inventory considerations 
include storage capacity, consumption characteristics and the consequences of stockouts at the cus- 
tomers. IRP’s first appeared in the literature in the 1970’s [2,12], however the bulk of the existing 
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research on problems of this type has appeared more recently (e.g. [1,3-11,13,14]). Practical optimization 
algorithms for realistically sized IRP’s have not been found and considerable ingenuity has been needed 
even to develop good heuristics. 

There are two different types of IRP’s — the Strategic IRP (SIRP) and the Tactical IRP (TIRP). The 
strategic version of the problem is motivated by the long lead times (months or even years) between the 
signing of purchase or lease agreements and the availability of vehicles for delivery operations. These 
long lead times make it impractical to build a fleet incrementally. SIRP focuses on estimating the 
minimum size (or cost) vehicle fleet required to supply inventory (e.g. industrial gases, fuel, consumer 
goods) from a central depot to spatially dispersed customers when only the probability distribution for 
the per unit time demand at each customer is known. TIRP, the related operational or tactical version of 
the problem, deals with routing an existing vehicle fleet to supply customers whose actual demands for 
replenishment are known or can be estimated. A major difference between the strategic and tactical 
versions of the problem is that in solving SIRP all possible realizations of TIRP must, at least implicitly, 
be considered simultaneously. 

The majority of work on IRP’s has focussed on TIRP. An exception is the SIRSA heuristic: developed 
by Larson [11] for a version of SIRP faced by the New York City Department of Environmental 
Protection in planning the acquisition of a barge fleet. SIRSA has a Clarke-Wright savings framework 
and estimates the minimum required fleet size by dividing the customers into a set of clusters. A basic 
assumption in SIRSA is that within any cluster, all replenishments are made on a single route visiting all 
customers in the cluster. When the solution developed for the New York City problem was examined, it 
was noticed that the Jamaica Bay location was visited several times a week, even though it needed to be 
visited only once every forty days. These frequent visits were the result of the inclusion in the cluster of 
other locations which needed to be visited frequently. Since a single route was used to service the entire 
cluster, the replenishment intervals were determined by the customer requiring the most frequent 
replenishments. The work reported here was motivated by this inefficiency and is based on'the use of 
period and phase of individual customer replenishment as additional decision variables. Period and 
phase of replenishment were first discussed in work on several versions of the Period Vehicle Routing 
Problem [3,10,12,13]. Introducing these variables enlarges the set of feasible replenishment schemes, so 
the resulting period/ phase heuristic is a generalization of SIRSA. 

The paper is organized as follows. Section 1 describes the SIRP model and reviews the SIRSA 
heuristic developed in [11]. Section 2 presents a simple example motivating the work in this) paper and 
defines the new decision variables — period and phase of customer replenishment. A petiod/ phase 
heuristic which exploits differences between the consumption and storage characteristics of individual 
customers is developed in Section 3. This heuristic is a generalization of SIRSA and uses the same 
savings framework. The relationship between our work on SIRP and earlier work on the PVRP is also 
discussed. Section 4 compares the performance of the period/ phase and SIRSA heuristics on a set of 
sixty test problems. The period / phase approach yields significant reductions in solution cost when the 
vehicle is large enough to replenish several customers in a single trip and/or when there is significant 
variation in the maximum inter-replenishment intervals of the customers. 


1. Problem statement and previous solution approach 


SIRP involves a single depot (with an inexhaustible supply of inventory) serving a set ¥ of N 
customers. All point-to-point travel distances d;, are deterministic and known. Customer i has inventory 
storage capacity C, and maximum acceptable probability of stockout between successive replenishments 
p;. If a customer stocks out, replenishment occurs on the next regularly scheduled visit to the customer. 


The time to load or unload q units of inventory at location i is given by a holding time function H,(q). 
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The daily consumption of inventory at any customer is a random variable with known mean m, and 

variance g,?. Demands on different days are assumed to be mutually independent, as are demands at 

different locations. The vehicle fleet will be composed of identical vehicles having travel speed V and 

carrying capacity W. The objective is to determine the minimum number of vehicles which should be 

acquired to service the customers such that 

(i) vehicle carrying capacity is not exceeded on any route, and 

(ii) the probability of any customer stocking out between successive replenishments is not greater than 
the customer’s maximum acceptable probability of stockout. 

The model above was appropriate for the problem faced in New York City. The assumptions used are 

similar to those in other inventory routing models and will be realistic for many applications. 


Previous approach for SIRP 
The SIRSA heuristic developed in [11] involves two steps: 


Step 1. Replace the original probabilistic demand SIRP with an equivalent deterministic demand 
problem DSIRP. 
Step 2. Solve DSIRP to find a fleet size estimate. 


Step 1: Define an equivalent deterministic demand problem DSIRP 

For customer i, define TMAX;, (the maximum inter-replenishment interval for customer i) to be the 
maximum possible time between successive replenishments (to full capacity C;) if the probability of 
stockout at customer i prior to the next replenishment is not to exceed p;. As described in [11] the 
cumulative demand at customer i over several days can be approximated by a normal random variable. 
Hence TMAX, must satisfy 


Pri Z> — | mp, (1) 
(TMAX,)o? 


where Z is a zero mean, unit variance normal random variable. Eq. (1) is readily solved for TMAX,. An 
effective daily consumption rate (or ‘worst-case’ daily demand) 4; > m, can then be calculated for each 
customer using 


yw; = C,/TMAX,. (2) 


t 


This effective daily consumption rate at customer i is the consumption rate which, if realized each day 
between successive complete replenishments TMAX, days apart, would exhaust the inventory of 
customer i at the instant the next replenishment occurs. 

DSIRP, an equivalent deterministic demand version of SIRP, is constructed by having customer i 
experience their effective daily inventory consumption of mu; every day. Using the effective daily 
consumption of each customer is a simple way of enforcing the chance constraints (1) without resorting 
to a complex mathematical model. 


Step 2: Solve the equivalent deterministic demand problem DSIRP 

The SIRSA fleet size estimate is based on a set of replenishment routes for DSIRP. These routes 
should, ideally, be optimal or near-optimal for all possible realizations of the tactical IRP the vehicle 
fleet may eventually meet. Since the number of such realizations is enormous, it is not practical to use a 
fully detailed model of the tactical problems in developing these routes. SIRSA generates a set of routes 
for DSIRP on the basis of a Simplified Tactical Model (referred to here as STM1) which is simple 
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enough to allow estimation of the fleet size but also retains many fundamental characteristics of the 

tactical IRP’s the vehicle fleet will eventually face. STM1 is based on the following four assumptions: 

(Al) The set of customers is divided into permanent disjoint clusters. 

(A2) Within a cluster, a single minimum length TSP route serves all customers. 

(A3) Customers are fully replenished whenever they are visited by a vehicle. 

(A4) A route is initiated only when necessary, i.e. when stockout is about to occur at a customer on the 
route or the accumulated demand for replenishment at the customers on the route is about to 
exceed vehicle carrying capacity. The second of these conditions forces route initiation because of 
assumption (A3). 

Assumption (A1) is included since any solution to an instance of the tactical IRP effectively divides 
the customers into clusters. The remaining assumptions define fleet operating policies which are 
consistent with maximizing the efficiency of vehicle use and estimating the minimum fleet size required. 
Assumption (A2) is an explicit statement of the standard operating policy in VRP’s and IRP’s. 
Assumptions (A3) and (A4) are similar to assumptions made in [6]. 

SIRSA solves DSIRP by developing a set of clusters (and associated routes) with minimum average 
vehicle requirement. For any given cluster Q, the average vehicle requirement p is defined to be 


t Total duration of the route serving Q 


es T Time between successive initiations of the route serving Q° 


Assumptions (A3) and (A4) result in route initiations at regular intervals of 


W 
T = mi in TMAX;,, ———}. 3 
min min j ca (3) 


ieN 


The first term in (3) ensures the maximum probability of stockout is not exceeded for any customer, 
while the second ensures that vehicle carrying capacity is not exceeded. The average vehicle requirement 
of a cluster can be interpreted as the number of vehicles needed, on average, to satisfy the replenishment 
requirements of the cluster. For any partition of “ into clusters the vehicle fleet size needed to service 
¥ is estimated by summing the average vehicle requirements of the individual clusters and rounding up 
the result to the next integer. 

The DSIRP clusters and associated fleet estimate are developed using a Clarke and Wright savings 
approach. The initial solution consists of a separate cluster and route for each customer. At each 
iteration, the saving (reduction in total average vehicle requirement) associated with each possible 
pairwise combination of clusters is calculated. The pair of clusters with the largest positive saving are 
merged. The process terminates when no combination of existing clusters will yield a positive saving. 

Adding the vehicle requirements of all clusters assumes that in actual daily operations each vehicle 
can be scheduled so as to be busy 100% of the time. Although this will not be true in practice, the fleet 
estimate obtained will be sufficient for actual scheduling problems due to the following sources of excess 
fleet capacity: 

(i) use of effective daily consumption at each customer means that in practice customers will not need 
to be visited as frequently as indicated in DSIRP; 

(ii) the savings algorithm is heuristic and will tend to overestimate fleet size; and 

(iii) on any particular day, the scheduler will be able to develop routes based on the current inventory 
levels of the customers. These routes will be more efficient for that particular day than the STM1 
routes on which the fleet estimate is based. 
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2. Period and phase of customer replenishment for DSIRP 


The SIRSA heuristic described above solves DSIRP by adopting simplified tactical model STM1. This 
section defines new decision variables — period and phase of customer replenishment — and describes a 
more general Simplified Tactical Model based on these variables. 


Motivation 


Consider the simple two-customer problem shown in Fig. 1. When STM1 is adopted, the feasible 
solutions are: 
(i) form separate clusters for customers 1 and 2 and use routes with inter-route intervals of, respectively, 
five and twenty time units, or 
(ii) form a single cluster and use the route 0-1-2-0 at intervals of five time units. 
These strategies result in the following average vehicle requirements: 
tor ttig — fo Hho 


PSIRSA (separate clusters) = Pcluster 1 + Petuster 2 = TMAX TMAX = 1.140 
1 2 


and 


to + bya + f29 Det 
Psresh combiner aus = Se aac GK ee ee: 
min{TMAX,, TMAX,} 5 


Now consider the following alternative service strategy: replenish customer 1 on days five, ten and 
fifteen, and on day twenty replenish both customers 1 and 2. The total travel time for all four routes is 
3(to, + 49) + (to + t12 + top) while the intervals between each successive pair of route initiations is five 
days. Repeated over time, this strategy satisfies the customer service and vehicle carrying capacity 
constraints and results in an average vehicle requirement of 


which is 16.2% lower than that for the best solution possible using STM1. The fleet estimate for the best 
STM1 solution is two vehicles while for the alternative solution it is one vehicle. 


Cust 1 
Cust 2 
Cust 0 
(Depot) 
TMAX, =5 toi = tio= 19 
TMAXp = 20 Travel Times ( tyg= tg, =2.0 
W>C + Co top = tog = 3.8 


Fig. 1. Two-customer inventory routing problem. 
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The inefficiency of the STM1 solutions is a result of assumption (A2) and the large difference in the 
maximum inter-replenishment intervals of the customers. This potential for inefficiency motivated the 
development of STM2 — a more general Simplified Tactical Model which allows differences between 
customer inter-replenishment intervals to be exploited. 


Period and phase of customer replenishment 


Before defining STM2, we introduce the concept of a routeset and define the new decision variables — 
period and phase of replenishment — used to characterize a routeset. 

For any cluster of customers, a routeset consists of a number of component routes arranged in a 
specific order. Each route is a minimum length tour which starts and ends at the depot and visits a subset 
of the customers in the cluster. The cardinality of a routeset is the number of routes it contains. A 
routeset is repeated over time to service the associated cluster. The alternative solution 


R, = (0-1-0, 0-1-0, 0-1-0, 0-1-2-0} 


suggested for the example in Fig. 1 is a routeset of cardinality | R, |= 4. 

The period of customer i (denoted by 7;) is the number of routes after a replenishment of customer i, 
up to and including the next route on which customer i is replenished. For example, in a cluster 
consisting of customers 1, 2 and 3 and serviced by the routeset, 


KR, = (0-1-3-0, 0-1-3-0, 0-1-2-3-0}, 


customers 1 and 3 have period one (they are serviced on every route) and customer 2 has period three (it 
is service on every third route). 

The phase of customer i (denoted by @¢;) is the number of routes between the beginning of a routeset 
and the first route replenishing customer i. In the routeset 


R, = {0-1-2-0, 0-1-3-0, 0-1-0, 0-1-0, 0-1-2-3-0}, 


customer | has a phase of zero and period of one, customer 2 has a phase of zero and period of four and 
customer 3 has a phase of one and period of three. 

Any routeset is completely characterized by its cardinality and the period and phase of each customer 
it services. 


STM2: A more general simplified tactical model 


The period/ phase algorithm described in Section 3 is based on the following Simplified Tactical 
Model, which we refer to as STM2. In STM2 each cluster is associated with a routeset rather than a 
single route. STM2 is defined by the following assumptions: 

(A5) The set of customers is divided into permanent disjoint clusters. 

(A6) Within a cluster, the customers are serviced by a single routeset. Each customer is visited by a 
non-empty subset of the component routes. 

(A7) Customers are fully replenished whenever they are visited by a vehicle. 

(A8) A route is initiated only when necessary, i.e. when stockout is about to occur at a customer on the 
route or the accumulated demand for replenishment at the customers on the route is about to 
exceed vehicle carrying capacity. 

Assumptions (A5)-(A8) are analogous to assumptions (A1)-(A4) in STM1. Since the single route 

allowed by (A2) in STM1 is a routeset of cardinality one, the feasible solutions for STM1 are a subset of 

the feasible solutions for STM2. 
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3. PPSA: The Period / Phase Savings Algorithm for SIRP 


The Period / Phase Savings Algorithm (PPSA) for SIRP has the same basic structure as SIRSA and 
involves the following two steps: 


Step 1. Replace the original probabilistic demand SIRP with an equivalent deterministic demand 
problem DSIRP using the equivalent daily demands yp,. 

Step 2. Solve DSIRP (using STM2) to find a period/phase solution with minimum average vehicle 
requirement. 


Step 1 is the same as for SIRSA. Step 2 employs a savings-type heuristic in which the basic computa- 
tional step involves finding, for any proposed cluster 2 of customers, a routeset with minimum average 
vehicle requirement. In SIRSA, Step 2 searches only among routesets of cardinality one. 


Minimum average vehicle requirement for a cluster 


The definition of average vehicle requirement given in Section 1 can now be extended to routesets. 
Given any cluster 2 and associated routeset R, let t; be the duration of route i and T, the interval 
between the initiations of routes i — 1 and i. T, is the interval between the initiation of the last route in a 
routeset and the initiation of the first route in the next repetition of the routeset. The average vehicle 
requirement for ¥, is defined to be 

|| |R| 
Pro= pay ae Be 8 (4) 
i=1t | i=] 
where the numerator is the total vehicle time required per repetition of the routeset and the denomina- 
tor is the time between successive initiations of the routeset. This definition is an extension of the 
definition in Section 1 since the feasible routes in STM1 are routesets of cardinality one. The average 
vehicle requirement for a routeset can be interpreted as the number of vehicles required, on average, to 
implement the routeset. 

We now describe procedure AvVehRegqt(), a two-stage heuristic which, for any specified Q, searches 
for a routeset with minimum average vehicle requirement. The first stage determines the cardinality of 
an initial incumbent routeset in which every route visits all customers in 2, while the second determines 
the period and phase of replenishment for each customer. 

Stage I: Determine cardinality of initial incumbent routeset. Given a cluster 2 containing n customers, 
this step determines the cardinality of the routeset , serving the cluster. Since the effort required in 
stage two increases with routeset cardinality, it is desirable to use as few routes as possible. Conversely, 
since each customer must be replenished at least once in the routeset, using too few routes may force 
some customers to be replenished more frequently than is necessary (recall the routesets of cardinality 
one in the example of Fig. 1). 

Assumption (A8) ensures that the initiation of a route is a result of either 
(i) the accumulated replenishment demand of customers on the route reaching the capacity of the 

vehicle, or 
(ii) the depletion of inventory at one or more of the customers on the route. 
From (3), a lower bound on the time between initiations of successive routes in the initial incumbent 
routeset is 


W 
LB = min{ min TMAX,, ; 5 
BE 2 o 


EN 
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Similarly, an upper bound on the maximum time between successive replenishment of any customer in 
the cluster is 


UB = min, max TMAX,, max me : (6) 
ieN ieQD ph; 
Since routes are initiated at least LB time units apart and the maximum time between successive 
replenishments of any customer is UB, the maximum number of routes between successive replenish- 
ments (and hence the maximum period) of any customer is bounded above by the ratio UB/LB. 
Employing a routeset of cardinality |§t,|=]UB/LB] will allow each customer to be replenished at least 
once in each repetition of the routeset without forcing unnecessarily frequent visits. 

Stage 2: Determine period and phase of replenishment for each customer in the cluster. We define an 
initial incumbent routeset for 2 consisting of | t,| routes, each visiting every customer in 2, i.e. t;=1 
and ¢,; = 0 for every customer i. This initial incumbent routeset is modified iteratively. At each iteration 
the customer/ period / phase modification (i*, 7*, #*) giving maximum reduction in average vehicle 
requirement is selected by a greedy search over all combinations (i, 7;, ¢;), 7;=2,..., | Hl and ¢,; = 0, 
1,...,7,-— 1, and for which 7;= 1 and ¢;=0 in the current incumbent routeset. The current incumbent 
routeset is updated by implementing (i*, 7*, ¢*). The iterative process continues until no modification 
reducing the average vehicle requirement is found, or until all customers have had their period and/or 
phase modified. 


Average vehicle requirement of a routeset: The route interval linear program 

A routeset is characterized by its cardinality and the period and phase of the customers it serves. 

Calculating the average vehicle requirement for a routeset involves the sequential calculation of: 

(i) the inter-route intervals T,, r= 1, 2,...,| Rol; 

(ii) the time intervals between successive replenishments of each customer (a function of the inter-route 
intervals and the customer periods and phases); 

(iii) the delivery quantity to each customer on each route and the total delivery quantity on each route; 

(iv) the total duration (travel time plus inventory transfer time) of each route; 

(v) the average vehicle requirement (using (4)). 

Steps (ii)-(v) are conceptually straightforward, although (iv) involves solving as many as |R,| 
travelling salesman problems. Calculating the inter-route intervals is not straightforward, as is illustrated 
by the routeset shown in Fig. 2. 

T,, the interval between route 3 (in the previous repetition of the routeset) and route | (in the next 
repetition of the routeset) is a function of the accumulated demand at customer 3, which is in turn a 
function of the intervals T, and T;. Similarly, customer 2 is visited only by route 2, so T, is a function of 
T, and T;. Since T, and T, are mutually dependent they must be calculated simultaneously. We now 
show how the inter-route intervals T, can be found by solving a linear program. 


Fig. 2. A simple three-customer, three-route routeset. 
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The inter-route intervals for the example in Fig. 2 are governed by the set of linear inequalities shown 
below. 
For route 1 (serving customers 1 and 3): 


T, < TMAX,, (maximum interval between replenishments of customer 1) 

T,+ T; + T, < TMAX,, (maximum interval between replenishments of customer 3) 

Tp, + (T,+T3+7,) 3 < W. (demand on route does not exceed vehicle capacity) 
For route 2 (serving customers 1 and 2): 

T, < TMAX,, (maximum interval between replenishments of customer 1) 

T,+7,+ T,<TMAX,, (maximum interval between replenishments of customer 2) 

Top, + (734+ 7, + Th) pe. < W. (demand on route does not exceed vehicle capacity) 
For route 3 (serving customer 1): 

T,; < TMAX,, (maximum interval between replenishments of customer 1) 

Tp, < W. (demand on route does not exceed vehicle capacity) 


Since the goal is efficient vehicle use L 7,, the interval between successive repetitions of the routeset, 
should be as large as possible. The problem of determining the route intervals can therefore be 
formulated as a Route Interval Linear Program (RILP) in which the T, are decision variables, the 
objective is to maximize & 7,, and the constraints are the linear inequalities above plus non-negativity of 
the 7,. 

In general, RILP is generated as follows. Let n be the number of customers in 2, | #, | the number of 
routes in the routeset and 2, Cf the set of customers served by route r (as determined by | #, | and the 
customer period and phase combinations being evaluated). For a maximization formulation the decision 


variables T,, r=1, 2,...,|9,| each have objective coefficient one. The constraints are of the form 
AT <b where T is the vector of decision variables. Constraint matrix A and right-hand side vector b 
have the structure described below. For route r, r= 1, 2,...,|Rgl: 


(a) There is a customer replenishment interval constraint for each i € 0,. Denote by (r, i) the row of A 
containing the constraint associated with route r and customer /. The value of a,, ;) ;, the coefficient 
of T; in this constraint, is 


1 for j=(L(r, i) +1), (LO, i) + 2),...r-1, 7, 
Q,.i),5 = where L(r, i) is the /ast route prior to r on which i was replenished, 
0 for all other j, 


ie. a,,,; is one if and only if customer 7 is not visited on routes j, j+1,...,r—1. 
The right-hand side value for this constraint is b,, ;, = TMAX,. 

(b) There is a vehicle capacity constraint. Denote (r, v) the row of A containing the vehicle capacity 
constraint for route r. Then a;,,,);, the coefficient of 7; in this constraint, is 


Qevj= Le My where M(r, j) ={iEO,: a, ,=1,1=1, 2,...,n}. 
kEM(r,j) 


M(r, j) can be interpreted as the subset of customers in 2, who are not replenished on routes j, 
j+1,...,7—1. For such customers, demand accumulated over the interval 7; must be met by the 
delivery on route r. 


The right-hand side value for this constraint is: b,,.,. = W. 


ryv) 
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If RILP is infeasible, the set of customer periods and phases being evaluated is also infeasible since it 
implies at least one negative element of T. 

RILP involves |" ,| decision variables and a maximum of |®,|-(7+ 1) constraints other than 
non-negativity constraints on T. There are at most n passes through Stage 2 and in each pass n+1—i 
customers are considered. For customer i, there are at most | %,|* different (7,, ¢,) combinations 
evaluated. Hence the number of instances of RILP solved in any call of AvVehReqt() is at most 
1) 9G |?-(n2 +n). 

A more formal statement of AvVehReqt() can now be given. 


Procedure AvVehReqt(Q) 

For a given incumbent routeset, let p, be the average vehicle requirement, U C2 the set of customers 
who have not had their period or phase altered, A, the largest possible reduction in p, if the period 
and/or phase of a single customer in U is altered and p,(i, 7;, 6;) the average utilization factor if the 
incumbent routeset is modified by changing the period and phase of customer i€U to 7; and ¢,; 
respectively. In addition, let 2, be the set of customers on route r, q, ; the inventory delivered by route r 
to customer j € 2,, Q, the total inventory delivered on route r, a, ; the time since the last replenishment 
of customer j € 2, L(r, j) the last route prior to r on which j was replenished and TSP(0,) the length 
of an optimal TSP tour starting and ending at the depot and visiting all customers in 2,. 


Stage 1. The cardinality of the initial incumbent routeset serving 2 is | ®,|=[UB/LB| where UB and 
LB are calculated using (5) and (6) above. Set U = 2. For each i € 2, set 7; = 1 and ; = 0. The 
initial incumbent routeset for Q has |#,| routes, each of which has O,=. For r=1, 


2,...,1 Kyl: 
T= T, = min{ min TMAX,, W/ > ui, 
ie ien 
9: =9,,=6;'T, ViEO, and Q=Q,= V q;. 


ieNn 


The average vehicle requirement for this initial incumbent routeset is 
pe (TSP(Q)/V + Ho(Q) + E H(a,))/T. 
feN 


Stage 2. If U#9: Set A, ~0. 
For each i € U: 
For 7; = 2, 3,...,|R pl and ¢;=0, 1, 2,...,7;-1: 
Calculate average vehicle requirement if period and phase of customer i are changed to 7; 
and @,. 
Form RILP. Solve for T. If RILP is infeasible, go to next i € U. 
For r= 1, 2,...,|R yl: Set OQ, =0. 
For each j © 2, calculate: 
a,j = ze T,., 
k=LUr,j)+1 
Gj Hj O, 5, 


Q, a QO; + 9, j° 
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The duration of route r is 


t, = TSP(2,) + H,(Q,) + yy Ai 4,4): 
keQ, 


The average vehicle requirement for this routeset is 


| Kt | Ft | 


pi, Ti» ¢;) = LG Di f,. 


r=1 r=1 


If p, — pi, 7, ;) > A,, set A, =p, — pi, 7;, 6), i* =i, T* = 7, and $* = ¢,. 

Stage 3. If A > 0, update the incumbent routeset by implementing the period / phase modification (i*, 
T*, o*). Set p, =p(i*, r*, $*) and U = U\ {i}. Go to Stage 2. 
If 4, <0, TERMINATE with the current incumbent routeset. The minimum average vehicle 
requirement for cluster 2 is p,. 


The Period / Phase Algorithm 


A formal statement of PPSA can now be given. A collection of routesets will be said to be a routeset 
cover if it services all customers without violating customer stockout or vehicle carrying capacity 
constraints. Let J be the number of clusters in a routeset cover and c the index for this routeset cover. 
Define p,; to be the average vehicle requirement associated with routeset j in routeset cover c. 


Step 1. Define DSIRP. For each customer i, calculate the maximum inter-replenishment interval TMAX;, 

and effective daily consumption yz; using Eqs. (1) and (2). 

Step 2. Solve DSIRP. 

Step 2a. Initialization of Savings Algorithm. The initial routeset cover consists of N routesets. Routeset j 

consists of a single one-stop route beginning at the depot, visiting customer j, and returning to 
the depot. Set J=N and c =0. Compute the average vehicle requirement p,, for routeset j: 
T(j) = Time between successive initiations of routeset j = min(W/u;, TMAX,}; 
t(j) = Total vehicle time required per repetition of routeset j 
= [do; + djgl/V + Ayu; TG) + Ayu; ‘T)); 
Po; = Average vehicle requirement for routeset j = ¢(j)/T()). 
Step 2b. Calculate savings. Consider pairwise each set of clusters j and k in the current routeset cover. 
For each pair of clusters, compute: 
(i) the minimum average vehicle requirement x, associated with retaining the two clusters as 
separate clusters: x, =p; + px; 
(ii) the minimum average vehicle requirement x, when the two clusters are combined. 
Let 2 be the set of customers if clusters j and k are merged. Use procedure AvVehRegt(2) 
to calculate x,; 
(iii) Compute the savings possible by combining clusters j and k: y,, = max{x, — x, 0}. 

Step 2c. Cluster modification. If at least one of the y,,’s computed in Step 2b is positive, select the largest 
yj, and revise the current routeset cover, replacing j and k by a single cluster involving the 
customers from clusters j and k. Eliminate from the previously computed p,,’s all those for 
which i =j or i=k. Reindex the routesets in the new routeset cover, increment c by one and 
decrement J by one. Return to Step 2b. 

If no y, is positive, TERMINATE: no further improvement is possible in the current routeset 
- cover. The fleet estimate is [Z/_1,,;l. 
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Table 1 

Sample problem data 

Location Holding Effective TMAX x-coord y-coord Setup Transfer 

index capacity consumption (days) (km) (km) time rate 
(units) (units /day) (hrs) (units /hr) 

0 oo 1 0 0.0 0.0 0.750 7000 

1 2000 1800 0.900 — 181.2 — 66.1 0.167 3500 

2 380 1690 4.447 — 180.5 60.4 0.167 3500 

3 680 2085 3.066 — 187.9 -11.9 0.167 3500 

4 1650 2030 1.230 — 143.5 152.7 0.167 3500 

5 3570 1515 0.424 — 199.3 — 14.6 0.167 3500 

6 2820 1560 0.553 — 36.2 25.3 0.167 3500 

7 630 1065 1.690 —174.1 16.4 0.167 3500 

8 1285 1880 1.463 — 120.6 162.3 0.167 3500: 

9 580 745 1.284 — 149.5 — 39.7 0.167 3500 


Illustrative example 

To illustrate the differences between the solutions generated by SIRSA and PPSA, the nine-customer 
problem in Table 1 was solved using both heuristics. Vehicle carrying capacity and travel speed are, 
respectively, 10000 units and 80 km/hr. Inventory transfer times consist of a setup time plus the time to 
transfer inventory at the indicated rate. The depot is location zero. The SIRSA solution is shown in 
Table 2. It consists of six clusters each served by a single route. The average vehicle requirement for this 
solution is 5.06, giving a vehicle fleet estimate of six vehicles. The PPSA solution consists of four clusters. 
The routeset details are shown in Table 3. The total average vehicle requirement for this solution is 4.73 
(a reduction of just over 6.5% cf. the SIRSA figure) and a vehicle fleet estimate of five vehicles (cf. six 
vehicles for SIRSA). 


Minimizing long-term average fleet cost 

PPSA has been developed in the context of minimizing the number of vehicles acquired. Minimizing 
the direct cost of fleet acquisition is a straightforward extension. Since the objective may involve 
minimizing some combination of acquisition and operating costs, we now indicate how PPSA can be 
extended to include approximate long-term operating costs per unit time. 

For a given cluster an associated routeset is a solution for STM2, which includes all major factors 
affecting the tactical IRP’s except actual customer inventory levels. The best routeset for a cluster can 
therefore be regarded as a set of ‘model’ routes which will be close to optimal for most tactical problems 
eventually encountered in the cluster. Even though the routes used in practice will vary to reflect actual 


Table 2 

SIRSA solution for sample problem — routeset details 

Cluster Route Customer sequence Length Duration Interval since Avge. vehicle 
(km) (hours) prior route (days) requirement 

1 1 0-1-5-0 447.2 77 0.424 2.267 

2 2 0-2-0 380.8 6.4 4.447 0.180 

3 3 0-7-3-9-0 408.8 74 1.284 0.721 

4 4 0-4-0 419.2 7.0 1.230 0.714 

5 5 0-6-0 88.3 2.7 0.553 0.608 

6 6 0-8-0 404.8 6.8 1.463 0.579 
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Table 3 
PPSA solution for sample problem — Routeset details 
Cluster Route Customer Length Duration Interval since Avge. vehicle 
index sequence h (km) (hours) prior route requirement 
(days) 
1 1 0-1-5-3-0 447.5 8.7 0.424 
2 0-5-0 399.7 6.6 0.424 
3 0-1-5-0 447.2 8.1 0.424 
4 0-5-0 399.7 6.6 0.424 
5 0-1-5-0 447.2 8.1 0.424 
6 0-5-0 399.7 6.6 0.424 
7 0-1-5-0 447.2 8.1 0.424 2.213 
2 1 0-7-9-0 390.8 6.6 1.284 
2 0-2-7-9-0 450.7 8.0 1.284 0.711 
3 1 0-4-0 419.2 7.0 1.230 0.714 
4 1 0-6-0 88.3 2.7 0.553 
2 0-6-0 88.3 2.7 0.553 
3 0-6-8-0 407.3 74 0.357 1.094 


customer inventories, the model routes can be used as indicated below to calculate approximations to the 
long-term operating costs incurred in supplying the customers in the cluster. 

For any specified cluster and routeset, the associated vehicle travel distance per unit time can be 
calculated from TSP(2,) and T,, the lengths and inter-route intervals for the component routes. Given 
travel cost data, the long-term travel cost per unit time for the routeset can be calculated and used as an 
approximation to the long-term travel cost per unit time for servicing the cluster. Similarly, the 7, give 
the lengths of the intervals between successive replenishments of each customer. Since the cumulative 
inventory consumption at customer i over several time periods is a normal random variable, the 
probability of stockout at customer i between any given pair of successive replenishments of i can be 
calculated. The expected number of stockouts per repetition of the routeset can be calculated for each 
customer. If a cost per stockout is specified for each customer, the long-term expected stockout cost per 
unit time can be calculated for the routeset. This can be used as an approximation for the actual 
long-term stockout cost per unit time for the cluster. 


Relationship between current work and the PVRP 


Period and phase of replenishment were first discussed in several papers (3,10,12,13] dealing with a 
vehicle routing problem for which the routing schemes used are similar in character to those generated 
in PPSA. This problem is referred to as the Period (or Periodic) Vehicle Routing Problem (PVRP) since 
the length of the planning horizon (or replenishment cycle) is a constant specified in the problem 
statement. 

There are significant differences between the work in the current paper and that in the papers cited 
above. The first is that the current paper is motivated by a different problem. While the four works cited 
above focus on the daily routing of an existing vehicle fleet over a specified planning horizon of several 
days, our work deals with fleet planning months or years in advance of actual daily operations. A second 
difference is that in [3,10,12,13] period and phase are measured in terms of numbers of time units. In our 
approach, period and phase are measured in terms of numbers of routes. A final difference is that in the 
model studied here, period and phase of customer replenishment are decision variables determined 
solely on the basis of vehicle carrying capacity, customer holding capacity and consumption characteris- 
tics, and the relative locations of the customers. In the PVRP papers the length of the planning period is 
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specified in the problem statement and restrictions are placed on the period and/or phase of 
replenishment. In [3], the period of replenishment for each customer is specified explicitly as a problem 
input and only phase of replenishment is a decision variable. In [13] a set of allowable k-day delivery day 
combinations for each customer is specified and the problem involves choosing one combination for each 
customer. Russell [12] studies a version of PVRP where the number of times a customer is to be 
replenished during the planning horizon is specified as an input. A similar restriction is used in [10] 
where it is assumed that each customer i is replenished exactly once in each interval of 7, days over a 
planning horizon of m days, where m/T; is assumed to be integer and m and T, are specified in the 
problem statement. These additional restrictions on period and/or phase of replenishment reduce the 
computation required to find a solution, but may also rule out potentially attractive solutions. 


4, Computational results 


The New York City data was unavailable. Testing of PPSA was carried out using sixty generated data 
sets. The test problems were designed to compare the performance of PPSA with that of SIRSA and to 
identify how any superiority of PPSA is related to two simple characteristics of the problem data: 

(i) RV, the ratio of the vehicle carrying capacity to the average customer holding capacity, and 

(ii) RT, the ratio of the largest maximum inter-replenishment interval to the smallest maximum 
inter-replenishment interval among all customers, i.e. 

RT= max {TMAX,}/ min {TMAX;}. 
i=1,2,...N i=1,2,...N 
These two ratios where chosen since: 

(i) PPSA is designed to exploit differences between the inter-replenishment intervals of customers. It is 
reasonable to expect that as these differences increase, PPSA should perform better relative to 
SIRSA. 

(ii) since there is usually a fixed cost associated with each stop a vehicle makes, a large cluster (i.e. one 
with many customers) is usually inefficient for small vehicles. For larger vehicles, larger clusters may 
be efficient. In larger clusters, it is more likely that the range of inter-replenishment intervals will be 
larger, with a corresponding improvement in the performance of PPSA relative to SIRSA. 

Two sets of base test problems were used. The first set (Uniform Data sets UD1—UDS) contained five 
base problems, each involving one hundred customers randomly located over a square distribution region 
of side 200 km centered on the depot. The second set (Clustered Data sets CD1 through CD5) contained 
five base problems, each having ten seed points randomly located in the distribution region. A cluster of 
ten customers randomly located in a square of side 50 km was centered on each seed point. For all base 
data sets, customer holding capacities were uniformly distributed between 5000 and 35000 gallons. The 
clustered problems were included since in actual problems customer locations may not be; uniformly 
scattered over the distribution region. An average TMAX value, TMAX,,, was chosen for each base 
data set. For base data sets one through five the TMAX,,-values chosen were, respectively, 2.5, 3.0, 3.5, 
4.0 and 4.5 days. 

Each base problem was solved for RV = 4 and 8 and RT = 5.0, 10.0 and 20.0, giving a total of sixty 
test runs. The TMAX-values of the customers were uniformly distributed in the range [2 - TMAX,,/(RT 
+ 1), 2-RT-TMAX,,/(RT + 1)] to give the required TMAX,,- and RT-values. Vehicle speed was 80 
km/hr and vehicle capacity was set equal to the product of RV and the average of the C;. The holding 
time functions used were 


H,(q) =0.75+1%10~°q for the depot, 
H,(q) =0.167+2%X107°q for the customers, 


where H(.) is in hours and g in gallons. 
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Table 4 
Computational results 
Base RT RV=4 RV=8 
data Avge. vehicle Imp. Fleet Time Avge. vehicle Imp Fleet Time 
requirement (%) redn. (seconds) requirement (%) redn. (seconds) 
SIRSA PPSA SIRSA PPSA’ SIRSA_ PPSA SIRSA PPSA 
UDI 5.0 8.0737 7.9928 10 1 912 10661 6.6879 6.4851 30 - 2054 12575 
10.0 9.6804 9.3265 3.7 - 662 18007 8.4389 7.8622 68 1 1941 64048 
20.0 11.766 11.608 13 - 927 11612 10.846 9.9177 86 1 1558 20353 
UD2 50 14.546 14.415 09 - 749 9064 11.710 11.760 0.2 - 1895 6451 
10.0 17.405 17.426 -01 - 921 10145 14.967 14.098 58 - 1847 38976 
20.0 21.968 20.939 47 1 794 21325 19.184 17.930 65 2 1271 84633 
UD3 5.0 8.5675 8.4512 14° - 1049 10885 7.1494 7.0164 19° = 3682 5820 
10.0 10.760 10.402 330 - 1378 15035 9.5424 9.0072 56 - 2049 29954 
20.0 14.273 13.513 53 1 1395 16886 13.358 12.042 99 1 2179 52851 
UD4 5.0 12.440 12.267 14° - 1164 8066 10.266 10.217 05 - 3228 6188 
10.0 15.394 15.047 23° - 831 11060 13.588 12.839 55 1 2159 19518 
20.0 18.217 17.738 2.6 1 1079 10892 16.850 15.246 95 1 2243 45271 
UD5 5.0 9.5080 9.4243 09° - 970 9985 7.8272 7.6807 19° ~ 2528 6968 
10.0 11.317 10.915 35 1 1110 14314 9.8594 9.2112 66 — 1768 19741 
20.0 13.915 13.109 58 - 948 19042 12.619 11.544 85 1 1624 31052 
CD1 5.0 74731 7.3935 110 - 769 23221 6.0842 5.9771 18 1 4269 15720 
10.0 9.0214 8.8316 2.1 1 865 14207 —s-_ 7.7891 7.4365 45 - 1560 26501 
20.0 11.128 10.951 16 1 1240 11762 =: 9.8536 )=—_- 9.4912 3.7 - 1525 60569 
CD2 5.0 14.209 13.967 17 1 645 16601 10.952 10.632 2.9 - 1318 17007 
10.0 17.260 16.943 18 1 844 28099 14.166 13.206 68 1 2247 78312 
20.0 21.755 20.651 5.2 1 1175 52247 = 19.848 18.098 88 1 2073 78964 
CD3 5.0 7.6714 7.6462 030 - 1005 17720 ~—- 6.2817 6.1317 2.4 - 1362 11508 
10.0 10.070 9.8189 25 1 1208 25902 =8.0517 ~—s_ 7.9925 0.7 1 2395 34771 
20.0 13.328 12.983 26 1 1527 22451 = 11.625 10.475 99 1 1986 98265 
CD4 5.0 11.693 11.471 19 - 779 20649 = 9.1786 =: 9.2294 -06 0 - 1633 7274 
10.0 14.377 14.019 25 - 799 15180 11.275 11.085 17) = 1838 14061 
20.0 17.405 17.082 19° - 887 14282 14.555 13.552 6.9 1 1796 42886 
CDS 5.0 8.6079 8.4893 14° - 969 14670 = 7.0140 ~— 6.7995 31 1 2089 10400 
10.0 10.600 10.483 11 - 784 16324 8.5472 8.1387 48 - 2443 30467 
20.0 12.655 12.303 28 - 1145 17397 10.996 10.104 81 - 1880 52721 


The results of the tests are shown in Table 4 and summarized in Table 5. The primary criterion on 
which solutions were compared is average vehicle requirement (which will usually be non-integer) rather 
than fleet size estimate (which is always integer). Average vehicle requirement is a better criterion for 
comparison since it is insensitive to the size of the fleets involved while fleet size (obtained by rounding 
up the average vehicle requirement) is not. For example, if the SIRSA average vehicle requirement is 


Table 5 

Average PPSA reduction in average vehicle requirement 

Problem RV=4 RV=8 

type RT=5 RT =10 RT = 20 RT=5 RT =10 RT = 20 


U 1.1 1.9 3.9 2.0 6.1 8.6 
Cc 
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99.5 even a small percentage reduction in average vehicle requirement will reduce the resulting fleet 
estimate by at least one vehicle. If however the SIRSA average vehicle requirement is 9.95 the fleet 
estimate will not change unless the PPSA average vehicle requirement is at least 9.5% lower than that 
for SIRSA. Computational times shown are in seconds on a Macintosh Ici desktop computer. 

PPSA produced reductions in average vehicle requirement in fifty-eight of sixty test problems. 
Reductions in fleet size were observed in twenty-six of sixty problems. As expected the performance of 
PPSA relative to SIRSA improves significantly as the variability in the maximum inter-replenishment 
intervals of the customers increases and/or vehicle capacity increases relative to average’ customer 
holding capacity. When RT = 20 and RV =8 the reduction in average vehicle requirement averages 
8.0% and fleet size reductions occur in eight of ten problems with these characteristics. When RT = 5 
and RV = 4 there is a 1% reduction in average vehicle requirement and fleet size is reduced in two of 
ten problems. 

It is worth noting that pursuing reductions in fleet size through the period/phase approach is 
worthwhile even in problems for which RT and/or RV are low. Reducing a fleet by even a single vehicle 
can result in substantial capital and operating cost savings (e.g. the vehicles involved in [11] wére barges 
costing approximately $5 million dollars each). 

PPSA performed somewhat better on the uniform problems than on the clustered problems, 
particularly for higher values of RT and RV. When the customers are clustered, the distance between 
the depot and the cluster will typically be much greater than the distances between customers within the 
cluster. If frequent visits are made to a customer with a small TMAX value, the additional cost'of visiting 
other customers in the cluster on each visit is relatively small, even if these other customers do not 
require frequent replenishment. For uniformly distributed customers, customer-to-customer: distances 
are larger relative to depot-to-customer distances and this additional cost is larger. As a result the 
‘one-route-per-cluster’ assumption of SIRSA is not as costly relative to PPSA when customer locations 
are clustered. 

The computational effort required by PPSA was considerably greater than for SIRSA. This is not a 
serious disadvantage since the potential benefits in term of reduced acquisition and operating costs for 
the fleet are much greater than the cost of using a desktop computer for the length of time required for 
PPSA. All but a handful of the test problems could have been solved overnight on a machine which is 
slow even by desktop computer standards. 
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